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This paper presents the analysis of axisymmetric helical magnetorotational instability (HMRI) in 
the inviscid limit, which is relevant for astrophysical conditions. The inductionless approximation 
defined by zero magnetic Prandtl number is adopted to distinguish the HMRI from the standard 
MRI in the cylindrical Taylor-Couette flow subject to a helical magnetic field. Using a Chebyshev 
collocation method convective and absolute instability thresholds are computed in terms of the 
Elsasser number for a fixed ratio of inner and outer radii A = 2 and various ratios of rotation rates 
and helicities of the magnetic field. It is found that the extension of self-sustained HMRI modes 
beyond the Rayleigh limit does not reach the astrophysically relevant Keplerian rotation profile not 
^v^j ■ only in the narrow- but also in the finite-gap approximation. The Keppler limit can be attained only 

by the convective HMRI mode provided that the boundaries are perfectly conducting. However, 
this mode requires not only a permanent external excitation to be observable but also has a long 
axial wave length, which is not compatible with limited thickness of astrophysical accretion disks. 
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>>. I. INTRODUCTION 

The magnetorotational instability (MRI) is a mechanism by which the magnetic field can destabilize a hydrody- 
namically stable flow of a conducting fluid without altering its velocity distribution. The MRI was first discovered 
^ . theoretically in cylindrical Taylor-Couette (TC) flow of perfectly conducting fluid subject to an axial magnetic field 
P, Three decades later, Balbus and Hawley @ suggested that the MRI may account for the fast formation of 
stars by driving turbulent transport of angular momentum in accretion disks. This proposition has triggered not only 
£>y numerous theoretical studies [4] but also several attempts to reproduce the MRI in the laboratory [a @|- A major 
f4 challenge to such experiments is posed by the parameter known as the magnetic Reynolds number Rm, which is re- 
quired to be at least ~ 10 for the MRI to set in. For a typical liquid metal, characterized by a small magnetic Prandtl 
number Pm ~ I0~ 5 — ICU 6 , this translates into a large hydrodynamic Reynolds number Re = Rm/Pm ~ 10 6 — I0 7 
£T} ■ 0|. At such large Reynolds numbers, the flow on which the MRI is expected to develop may become turbulent due 
J> ' to purely hydrodynamic instabilities [f|. 

A way to circumvent this problem was proposed by Hollerbach and Riidiger [8], who suggested that a 
magnetorotational-type instability can take place in cylindrical TC flow at Re ~ I0 3 when the imposed magnetic 
field is helical rather than purely axial as for the standard MRI (SMRI). An instability resembling this new type 
helical MRI (HMRI) was shortly thereafter observed in the PROMISE experiment [9KL1|. These observations were 
disputed by Liu et al. [12], who found no such instability in their inviscid theoretical analysis of finite length cylin- 
ders with insulating end caps. Subsequen tly, the observed phenomenon was conjectured to be a transient growth 
rather than a self-sustained instability (f3l - tl5( . This proposition was, in turn, opposed by Priede and Gerbeth [llij . 
who showed that there is not only a convective but also an absolute HMRI threshold. Thus, a self-sustained HMRI 
can experimentally be observed in a system of sufficiently large axial extension. However, the comparison with the 
experimental results revealed that the HMRI has been observed slightly beyond the range of its absolute instability, 
where it is expected according to ideal TC flow model. This discrepancy implies a deviation of the real base flow from 
the TC one, which may be caused by the end effects. One such end effect disturbing the base flow in the original 
PROMISE experiment was due to the copper end cap, which allowed radial electric current to connect over the liquid 
gap between the cylinders [lj}- This flaw was corrected in the modified PROMISE experiment, where an insulating 
end cap was used, which, in order to reduce the Ekman pumping, was split into two separately rotating rings [iH , [l^ | . 
Although the instability appeared much sharper in the modified setup than in the original one, its nature may still be 
questionable. Uncertainty is due to the virtually unknown hydrodynamic stability limit of the real base flow, whose 
deviation from ideal TC flow was significant and practically unavoidable in the experiment. Without knowing the 
actual hydrodynamic limit it is practically impossible to distinguish the HMRI from a magnetically modified Taylor 
vortex flow [I7fl ■ 
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Figure 1: (Color online) Sketch of the formulation of the problem. 

Regardless of the experimental reproducibility, the HMRI is also of a questionable astrophysical relevance. Using 
a WKB analysis Liu et al. [HJ showed that however the HMRI is able to destabilize centrifugally stable velocity 
distributions, it does not reach up to the astrophysically relevant Keplerian rotation profile. This claim was doubted 
by Riidiger and Hollerbach |20|, who pointed out that according to the numerical results for the TC flow in a 
finite- width gap [2lJ, the HMRI can apparently reach the Keplerian rotation profile provided that at least one of the 
boundaries is sufficiently conducting. However, it is important to note that this conclusion concerns only the convective 
HMRI, which is not generally self-sustained and requires an external excitation to be effective. No destabilizing effect 
due to conducting boundaries was observed for the absolute instability threshold [16]. Although the extension of the 
absolute HMRI beyond the Rayleigh line was found to increase with the strength and helicity of the magnetic field, 
the previous study was unable to conclude whether the Keplerian velocity profile is attainable. The problem was the 
excessive numerical resolution required for the thin boundary layers developing in strong magnetic field. However, the 
HMRI does not appear to be related with the boundary layers, whose main function is to satisfy the no-slip boundary 
condition imposed by the viscosity. These redundant boundary layers can be eliminated by ignoring the viscosity, 
which appears to be insignificant for astrophysical conditions. In this paper, fluid is considered not only inviscid but 
also highly electrically resistive. It means that the magnetic Prandtl number is assumed to be zero regardless of the 
viscosity, which corresponds to the inductionless approximation [2l| . 

The aim of the present work is to investigate numerically whether the astrophysically relevant Keplerian rotation 
profile can be attained by inviscid HMRI in cylindrical TC flow when a finite-width annulus is considered. The 
obtained results show that only convective HMRI mode can reach the Kepler limit provided that the boundaries are 
perfectly conducting. However, this mode requires not only a permanent external excitation to be observable but also 
has a long wave length, which is incompatible with limited thickness of accretion disks. The absolute HMRI as well 
as the convective one at insulating boundaries are found to obey the Liu limit also for a finite-width annulus. This is 
because the maximum extension of the HMRI beyond the Rayleigh line is attained when the magnetic field is nearly 
azimuthal, which results in a short-wave instability as in the local WKB approximation. 

The paper is organized as follows. In Sec. [XT] the problem is formulated. The local WKB-type solution for a 
narrow gap is revisited in Sec. IIII1 Section IIVI presents numerical results concerning the convective and absolute 
HMRI thresholds for both insulating and perfectly conducting boundaries. The paper is concluded with a summary 
of results in Sec. [V] 

II. PROBLEM FORMULATION 

Consider an incompressible inviscid fluid of electrical conductivity a filling the gap between two infinite concentric 
cylinders with the inner radius Ri and the outer radius R a rotating with the angular velocities Qi and f2 Q , respectively, 
in the presence of an externally imposed steady magnetic field B = B^e.^ J rB z e. z with axial and azimuthal components 
B z = Bq and B^ = (3BoRi/r in the cylindrical coordinates (r, </>, z), where /3 is a dimensionless parameter characterizing 
the geometrical helicity of the magnetic field (see Fig. 1). The induced magnetic field is assumed to be negligible 
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relative to the imposed one, which corresponds to the inductionless approximation. The velocity of inviscid fluid flow 
v is governed by Euler equation with electromagnetic body force 

<9v i 

- + (v.V)v = p- 1 (-Vp+jxB ), (1) 

where the induced current follows from Ohm's law for moving medium 

j = a(E + vxB ). (2) 

In addition, we assume the characteristic time of velocity variation to be much longer than the magnetic diffusion time 
to 3> T m = noaL 2 , which leads to the quasi-stationary approximation, according to which VxE = and E = — V$, 
where <I> is the electrostatic potential [22j. Mass and charge conservation imply V • v = V • j = 0. 

As long as the viscosity is not exactly zero, the problem admits a base state with a purely azimuthal velocity 
distribution vq(t) — e0Do(r), where the velocity profile 



v {r) = r 



Rl-Rj ' rR~ 2 -R~ 2 



is independent of the viscosity and thus holds also in the inviscid limit. 

Note that the base flow is not affected by the magnetic field because the latter gives rise only to the electrostatic 
potential &o(r) = Bq J vo(r)dr, whose gradient compensates the induced electric field so that there is no current in 
the base state (j = 0). However, current may be induced in the perturbed state 

ri}(-«»-{2:«}M + { r.«}(-" 

where Vi, pi, j 1; and $i are small-amplitude perturbations for which Eqs. (TTJ) and @ after linearization take the 
form 

^ + (vi • V)v + (v • V)vi = p- 1 (-Vpi + Ji x B ) 

j x = o-(-V$ x +Vi x B ). (3) 

In the following, we focus on the axisymmetric perturbations which are typically much more unstable than the 
nonaxisymmetric ones [23j. The nonaxisymmetric HMRI modes recently found by Hollerbach, Teeluck and Riidiger 
[24j in purely azimuthal magnetic field lay outside the scope of this study. For axisymmetric perturbations, the 
solenoidality constraints are satisfied by meridional stream functions for fluid flow and electric current as 

v = ve^ + V x ^e^, 
j = je + V x he$. 

Note that h is the azimuthal component of the induced magnetic field, which is used subsequently instead of $ for 
the description of the induced current. Thus, we effectively retain the azimuthal component of the induction equation 
to describe the meridional components of the induced current, while the azimuthal current is explicitly related to 
the radial velocity. The use of the electrostatic potential which provides an alternative mathematical formulation 
for the induced currents in the inductionless approximation, would result in slightly more complicated governing 
equations. In addition, for numerical purposes, we introduce an auxiliary variable, vorticity 

ijj = V x v = cje^ + V x ve^. 

Perturbation is sought in the normal mode form 

Wi.Vi, hi} (r, t) = {«, u>, h} (r) x e^ t+lkz , (4) 

where 7 is, in general, a complex growth rate and k is the axial wave number, which is real for the conventional 
stability analysis and complex when the absolute instability is considered. Henceforth, we proceed to dimensionless 
variables by using Ri, Ri^li, Bq, and aB^Ri^li as the length, time, velocity, magnetic field, and current scales, 

respectively. The nondimensionalized governing equations read as 

71) = ikr^ 1 (r 2 rt)'ip + Aikh, (5) 

7cj = 2ikQv~ Aik(iktp + 2/3r' 2 h), (6) 

= D k ^ + uj, (7) 

= D k h + ik(v-2Pr~ 2 T]j), (8) 
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where D k f = r" 1 (rf) - (r~ 2 + k 2 )f and the prime stands for d/dr; A = aB 2 /pili is the Elsasser number, also 
referred to as the magnetohydrodynamic interaction parameter, which will be the main dimensionless parameter for 
the subsequent analysis; 

is the dimensionless angular velocity of the base flow defined in terms of A = R /Ri and fj, = fi /flj. Note that Eqs. 
©-© are invariant upon the transformations {SI, k, h} —> — {fi, k, h} and {/3, k, v} —> — {/?, k, v}, which means that 
reversing the direction of rotation Q or that of helicity /3 merely inverts the direction in which perturbations propagate 
along the axis. Similarly, Eqs. ©-([5]) are invariant upon the transformation [k, 7) — > (—k, 7*) applied together with 
the complex conjugate operation, which implies that for a fixed wave number the reversion of SI or /3 inverts only the 
sign of frequency u) — Impy]. Owing to these symmetries it suffices to consider only positive 

In inviscid approximation, the flow perturbation on the inner and outer cylinders at r = 1 and r = A, respectively, 
satisfy the impermeability condition ip = 0. Boundary conditions for h on insulating and perfectly conducting cylin- 
ders are h = and (rh)' = 0, respectively. There are no boundary conditions imposed on the azimuthal velocity 
perturbation, which is defined in terms of the two previous variables by Eq. 

Equations ([5])-([5} for perturbation amplitudes were solved in the same way as in Refs. [ill [2lj using a spectral 
collocation method based on the Chebyshev-Lobatto grid with a typical number of internal points N = 32, which 
ensured the accuracy of about five digits. 



III. NARROW-GAP APPROXIMATION 



Before undertaking numerical solution of the whole problem, it is instructive to revisit the local WKB-type solution 
based on the narrow-gap approximation defined by <5 = A — 1 < 1, which is the dimensionless gap width between the 
cylinders. Representing the cylindrical radius as r = 1 + 5s, where s is a local radial coordinate, the angular base 
flow velocity © in Eq. (jX3J) at the leading order in 8 becomes SI ps 1. The dimensionless radial gradient term of the 
angular momentum in Eq. ([S]) can be written for arbitrary gap width as 

r-VS!)' = 2^-5— |=2(l + Ro), (10) 

which is constant not only for narrow but also finite-width gap. The latter relation can be used to define the Rossby 
number for arbitrary gap width as 

Ro= F7 ^ T , (11) 

which for a narrow gap reduces to Ro s» (/i — l)/(28) « (ln-v/f2)'. Taking d as the length scale, which implies the 
substitutions {k, j3} —> {k, /3}/5, {-0, h} —> {ip, h}5 2 and v — > vS, Eqs. ©-© reduce at the leading order to 

76 = 2ik(l + Ro)ip + Aikh, (12) 

7w = 2ikv- Mk(ik4> + 2l3r- 2 h), (13) 

= {jf -k 2 ^ + uo, (14) 

= h" -~k 2 h + ik{v-2l3r- 2 rjj), (15) 

where the prime now stands for d/ds. 

A locally periodic solution can be sought in the form {v, ip,uj,h}(s) = {vo, ipoyiuo, ho} sin(Zs) with constant amplitudes 

and the radial wave number /. Substituting this into Eqs. (|12l) - (|15[) and eliminating the amplitudes yields the 
dispersion relation 

4 (ik(l + Ro) + fiku 2 ) (ik + f3An 2 ) - (7 + An 2 ) ((7 + Ak 2 )K 2 + 4/3 2 Ak 2 ) = 0, (16) 

where K 2 — k 2 + I 2 and k — k/K. 

Although the quadratic equation ([TO)) is easy solvable for the growth rate 7, it is instructive first to consider 
asymptotic solutions for small and large interaction parametersA. In the nonmagnetic case A = 0, the solution is 
simple: 

70 = ±2 Wl + Ro, (17) 
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which describes exponentially growing Taylor vortices due to the centrifugal instability when Ro < — 1, and constant- 
amplitude inertial waves when Ro > — 1. For strong interaction at A 1, Eq. ([To]) , reducing to 



4(/3Aac 2 ) 2 - ( 7 + An 2 ) (( 7 + Ak 2 )K 2 + 4(3 2 An 2 ) = 0, 



yields 



7 = -Ak 2 1 + 2(0/ K f 1 ± y/1 + (p/K) 



<0, 



(18) 



which shows that all perturbations are magnetically damped. For weak interaction at A <C 1, the growth rate may be 
sought as a perturbation of the nonmagnetic solution 7 rs 70 + A71 , where 70 is defined by Eq. (|17|) . This expression 
substituted into Eq. (fTo| yields 



7i 



-K' 



k A + 2(/3k) 2 ± fink 



Ro 



Vi + Ro 



According to the above expression, there may be a finite range of helicities bounded by 

2 



K 

T 



Ro 



VT + Ro 



± 



(2 + Ro)- 
1 + Ro 



(19) 



(20) 



in which inertial waves turn unstable (71 > 0) however small the interaction parameter A. It means that the inviscid 
HMRI, similarly to its standard counterpart, can be triggered by an arbitrary weak magnetic field provided that its 
helicity lies in the range defined above. For such an instability range to exist, /3± has to be real, which, in turn, 
requires the term under the square root to be nonnegative and yields 



±Ro > 2V2± 2 = ±Ro^. 



(21) 



Taking the plus and minus signs in this expression yields Ro > 2 + 2V2 = RoJ and Ro < 2 — 2y/2 — Ro^, respectively, 
which are originally due to Liu et al. [12]. 

The damping of all perturbations at A 1 following from Eq. (118p implies that instability, if any, is possible only 
at sufficiently small A. The marginal interaction parameter A c , at which perturbations become neutrally stable, is 
defined by Re [7] = 0. Substituting 7 = ioj, where uj is the frequency of neutrally stable perturbations, into Eq. |T| 
and taking the imaginary part, we obtain 

2(2 + Ro)k 



K//3 + 2p/K' 

This expression substituted into the real part of Eq. (fT5|) results in 



A, 



V. 2 /. 2 - 4(1 + Ro) = 2^ ( K ^; /K ) - (1 +"^)- 



(22) 



Again, for instability to exist, A c has to be real, which requires the term under the square root to be non-negative. 
This is equivalent to 



> (K/P + 2(i/K) 2 = f(K/P) > f(V2) 



(2 + Ro) 2 
1 + Ro 

which results in the same constraint as that for A <C 1 given by Eq. (|21l) . When the condition above is satisfied, 
instability is confined to the finite wave number band defined by the limiting values, 

(23) 

^ 1 yit no v 1 -r / 

at which A c in Eq. ([22]) turns to zero. 

These are the basic characteristics of the HMRI in the narrow-gap approximation. It is important to note that 
although this instability can affect some centrifugally stable velocity distributions with Ro > —1, it does not extend 
up to the astrophysically relevant Keplerian profile Q ~ r -3 / 2 , whose Rossby number (ITU) . 

3 




Ro K 



> Ro] 



-0.828, 



(24) 



lies outside the Liu range ([2"T]) . The aim of the following section is to investigate numerically whether this constraint 
can be overcome when a finite-width gap is considered. 
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Figure 2: (Color online) The marginal interaction parameter A against the wave number k for the helicities (a) /3 — 2 and (b) 
f3 — 8 at different /i; A = fi/k at the Rayleigh line /i = 0.25. 

IV. NUMERICAL RESULTS FOR A FINITE- WIDTH GAP 

The numerical results presented in this section are for the gap width between the cylinders equal to the radius of the 
inner cylinder, which corresponds to A = R /Ri = 2. The corresponding Rayleigh line is hr — o /^i = A -2 = 0.25. 
Above this critical ratio of rotation rates of the inner and outer cylinders the specific angular momentum r 2 Cl turns 
radially outward increasing and thus the flow becomes hydrodynamically stable with respect to axisymmetric inviscid 
perturbations. 

A. Convective instability 

In this section, real wave numbers are considered (fcj = 0) as in the standard linear stability analysis, which corre- 
sponds to the so-called convective instability. At this threshold, flow becomes able to amplify certain perturbations, 
which are not necessarily self-sustained and thus may require a permanent external excitation to be observable. For 
self-sustained perturbations, the absolute instability is required, which is considered in the next section. 

The marginal interaction parameters below which the temporal growth rate j r becomes positive, are shown in 
Figs. [2] and [3] against the wave number of the corresponding mode for insulating and perfectly conducting boundaries, 
respectively. For the centrifugally marginal state defined by fi — the neutral stability boundary turns out to be 
particularly simple: A = f3/k. The origin of this simple relation is not obvious. Slightly beyond the Rayleigh line 
(fi > 0.25), the short-wave unstable modes disappear. With the increase in fi, the range of unstable modes rapidly 
shrinks toward the longer waves. As seen in Fig. [5Ja), for sufficiently low helicities /3, there is a critical \i at which the 
long-wave instability mode disappears altogether. For higher /3, the instability disappears in a more complex way. As 
seen in Fig. HJb) for j3 = 8, there is an intermediate-wave mode, which outlasts the long-wave one. This longer lasting 
mode affects only a limited range of wave numbers, which also quickly shrinks from both ends with the increase in 
(i. As a result, the range of unstable wave numbers disappears altogether by shrinking to a point at some critical 
(J,, which slightly exceeds that for the disappearance of the long-wave mode discussed above. As seen in Fig. [H the 
critical wave number at which this mode disappears increases with the helicity j3. Thus, this intermediate-wave mode 
resembles the narrow-gap instability considered in the previous section. 

For perfectly conducting boundaries and sufficiently low helicities, the marginal interaction parameter, which is 
shown in Fig. [S^a) against the wave number for /3 = 2, demonstrates a similar variation to that for insulating 
boundaries seen in Fig. [Ha). In contrast to the insulating boundaries, with the increase in //, the original long- 
wave instability mode quickly disappears. What is left behind resembles an intermediate- wave mode for insulating 
boundaries at high helicities except that the remaining mode extends up to the long- wave limit k = 0. In contrast to 
the original long- wave mode at /j, — the marginal interaction parameter for this mode tends to a finite value rather 
than growing as ~ k^ 1 for k — > 0. The main difference from the insulating case is that instability does not vanish 
at a finite /i. With the increase in instability just shifts to longer waves with the marginal interaction parameter 
approaching zero, which, as shown in the previous section, implies a similar decrease also in the growth rate. The 
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Figure 3: (Color online) The marginal interaction parameter A against the wave number k for the helicities (a) /3 — 2 and (b) 
f3 — 16 at different /i; A = f3/k at the Rayleigh line /i = 0.25. 
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Figure 4: (Color online) Extension of the HMRI beyond the Rayleigh line (fi = 0.25) for (a) insulating and (b) perfectly 
conducting boundaries at various helicities /3. 



variation of the marginal interaction parameter for perfectly conducting boundaries at higher helicities is shown in Fig. 
[3]^b) for j3 — 16. Slightly beyond the Rayleigh line, instability is due to an intermediate-wave mode, which vanishes 
at a finite /i similarly to that for the insulating boundaries seen in Fig. [Hb). What is left behind is a long- wave 
mode, which persists at large (x by shifting to ever longer waves as in the case of low helicities considered above. The 
above results are summarized in Fig. |4j which shows the extension of instability beyond the Rayleigh line in terms of 
the marginal wave numbers at which the instability vanishes plotted against the ratio of rotation rates [i for various 
helicities (3. 



B. Absolute instability 



The convective instability considered in the previous section is not necessarily self-sustained. In spatially unbounded 
systems, small amplitude perturbations become self-sustained above the absolute instability threshold. In addition to 
the positive temporal growth rate supposed by the convective instability, the absolute instability requires zero group 
velocity for the critical perturbation. This additional constraint can be satisfied by a nonzero imaginary part of the 
wave number fc, ^ 0. From a physical as well as computational point of view, it is advantageous to consider the 
absolute instability as a limiting case of the global instability. The latter supposes that two oppositely propagating 
modes with the same frequency can be mutually coupled by the reflections from remote end walls. Two such marginal 




Figure 5: (Color online) Interaction parameter A (a) and the frequency w (b) for marginally stable modes (7,. = 0) versus the 
real part of the wave number k r for [i — 0.27 and f3 = 5 at various imaginary parts of the wave number fe. 



modes having the same imaginary but, generally, different real parts of the wave number can form a neutrally stable 
global mode with a finite wave packet length fitting in a system of sufficiently large extension |25l . |26| . The global 
instability turns into the absolute one when two such oppositely propagating modes merge so that also the real parts 
of the wave number coincide, which results in a wave packet of infinite extension. This approach will be pursued in 
the following to locate the absolute instability. 

The first step is to consider the neutral stability curves for fcj / which are plotted in Fig. [5] for insulating 
boundaries with a rotation rate ratio slightly beyond the Rayleigh line (/x = 0.27) and a moderate helicity (3 = 5. 
As discussed in the previous section, the neutral stability curve for ki = 0, which defines the convective instability 
threshold, consists of an intermediate- and a long-wave branch. Transition between these branches shows up as a 
break on the neutral stability curve and as a jump in the associated frequency, which are seen in Figs. [SJa) and 
EJb), respectively. For ki > representing the main interest here, both branches merge together and neutral stability 
curves start to form closed loops, which tighten as ki is increased. A more relevant information follows from Fig. [5J 
where the previous two quantities are plotted against each other. It is important to notice that the curves in Fig. 
|6^b) form not only closed loops but also intersect themselves in a certain range of ki. For a fixed ki, the frequencies 
and the interaction parameters of two modes coincide at the intersection. As discussed above, two such waves could 
sustain each other by the reflections from the end walls and thus form a marginal global mode provided that they 
propagate in the opposite directions. The direction of propagation can be determined from the variation of the 
interaction parameter with ki. Namely, two modes propagate in the opposite directions if the interaction parameter 
of one intersecting branch increases while the other decreases with a small variation in ki (26| . Since this is the case 
seen in Fig. |H{b), the intersection point defines a marginal interaction parameter for a neutrally stable global mode 
consisting of two oppositely traveling waves with different real parts of the wave number. At ki = 3.80, the loop below 
the intersection tightens together forming a cusp. This critical point at which both modes merge together represents 
the threshold of the absolute instability. Besides the lower threshold there is also an upper one, which occurs at 
ki = 3.28 when another cusp forms at the top of the loop. This upper cusp is formed by the second intersection which 
is hardly noticeable at the upper tip of the loop for ki = 3 in Fig. [f^b) . 

The thresholds of the absolute instability for insulating boundaries at various helicities are plotted in Fig. [7] against 
the ratio of rotation rates \i. With increasing the rotation rate or reducing the strength of the magnetic field, the 
HMRI sets in below the upper critical value of the interaction parameter and vanishes below the lower one, which 
are shown in Fig. EJa). Such a double threshold is a characteristic feature of the HMRI that distinguishes it from 
a magnetically modified Taylor vortex flow, which exists below the Rayleigh line and has only the upper threshold. 
The critical phase velocity and the real and imaginary parts of wave number at the instability thresholds are shown 
in Figs. EIb)-[7{d), respectively. Although the increase of the helicity j3 makes the instability extend farther beyond 
the Rayleigh line, the extension is seen in Fig. [9] to saturate at large /3. 

For the absolute instability, in contrast to the convective one considered in the previous section, there is no significant 
difference between insulating and perfectly conducting boundaries. For perfectly conducting boundaries, only the 
intermediate but not the long- wave modes can turn absolutely unstable. As seen in Fig. [SJa), the neutrally stable 
modes with ki > are limited to finite wave number and interaction parameter ranges. The marginal interaction 
parameter plotted in Fig. \Eib) against the frequency of neutrally stable modes forms self-intersecting loops as in 




Figure 7: (Color online) The critical interaction parameter A c (a), the phase velocity (b), the real (c) and imaginary (d) parts 
of the critical wave number of the absolute instability versus fi for insulating boundaries at various helicities j3. 
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the case of insulating boundaries shown in Fig. [U[b). Thus, the threshold of the absolute instability for perfectly 
conducting boundaries turns out to be close to that for insulating boundaries shown in Fig. [7J 

The proximity of the absolute instability thresholds for insulating and perfectly conducting boundaries is demon- 
strated in Fig. [5J which shows the extension of the HMRI beyond the Rayleigh line depending on the helicity /?. Note 
that the extension of the convective instability for perfectly conducting boundaries is not shown because it appears 
unlimited in Fig. 01b). At high /?, the extension of the convective instability for insulating boundaries as well as that 
of the absolute instability for both insulating and perfectly conducting cylinders are seen to tend to a certain limit of 
/i. The main question is whether this limit includes the Keplerian velocity profile Q ~ r -3 / 2 . A simple approximation 
of the Keplerian profile with the TC one leads to /ia = A~ 3 / 2 « 0.35 for A = 2 under consideration here. As seen in 
Fig. |2 the HMRI certainly persists beyond this limit. On the other hand, the lower Liu limit (|2 1 1) substituted into 
Eq. IQIJ yields 

Mi = l + (l-A- 2 )Ro^^--^=, (25) 

which appears to be the asymptote approached by /x in Fig. [HI It means that the HMRI obeys the Liu limit not only 
for a thin but also for a finite-width gap. But this obviously contradicts the above conjecture that the HMRI can 
affect the Keplerian profile, whose Rossby number according to Eq. (|2~3| lies beyond the Liu limit and yields 

UK = 1 + (1 - A- 2 )Ro x = — . (26) 

16 

Indeed, this ratio of rotation rates is seen in Fig. [H] to lie significantly above that for the Liu limit. It is important 
to note that /za is an integral criterion based on the difference of angular velocity over the gap width. The velocity 
difference for the TC profile coinciding with that for the Keplerian profile means that only the average velocity 
gradients of both profiles are the same. Since the velocity profiles are different, the TC profile will be in some 
parts of the gap shallower and in some other parts steeper than that of the Keplerian profile [20] . Therefore, simple 
approximation of the Keplerian profile by the TC one based on the total velocity variation may be misleading. This 
becomes obvious by considering the radial gradient term of the specific angular momentum given by Eq. (j 10[) . which 
determines the centrifugal stability. This is also the key parameter for the HMRI, whose critical value is modified 
by helical magnetic field as shown by the narrow-gap approximation. Expression (|10p defines the local slope of the 
specific angular momentum profile, which is constant for TC flow and can be expressed in terms of either fi or Ro. 
Then it is easy to see that the slope of TC flow profile corresponding to /ia is indeed considerably lower than the 
local slope of the Keplerian profile defined by Eq. (|26l) . 



V. SUMMARY AND CONCLUSIONS 



This paper presented the analysis of the HMRI in the astrophysically relevant limit of a vanishing viscosity. In 
order to focus on the HMRI, the analysis was carried out in the inductionless approximation defined by zero magnetic 
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Figure 9: (Color online) The extension of the convective and absolute HMRI beyond the Rayleigh line versus the helicity ft for 
insulating and perfectly conducting boundaries. 

Prandtl number Pm = 0, which excludes the standard MRI. These assumptions resulted in the problem defined by 
a single dimensionless parameter, the Elsasser number, also referred to as the MHD interaction parameter, which 
characterizes the ratio of electromagnetic and inertial forces. Using a Chebyshev collocation method convective and 
absolute instability thresholds were numerically calculated with respect to axisymmetric perturbations of cylindrical 
TC flow subject to a helical magnetic field. Without the magnetic field, the flow becomes centrifugally unstable when 
the specific angular momentum of fluid at the inner cylinder exceeds that at the outer one. According to the previous 
WKB analysis, a strong magnetic field has a stabilizing effect, but a weak field of certain helicity can destabilize 
some centrifugally stable velocity distributions lying beyond the so-called Rayleigh line. The latter defines the limit 
of the centrifugal instability in terms of the critical ratio of rotation rates for a given gap width between the cylinders. 
However, in the WKB approximation, the extension of the instability beyond the Rayleigh line is constrained by the 
so-called Liu limit and does not reach up to the astrophysically relevant Keplerian velocity profile. The obtained 
numerical results show that the Liu limit is obeyed also by the nonlocal solution for a finite width gap between the 
cylinders except for perfectly conducting boundaries. Maximal extension of the absolute HMRI as well as that of the 
convective HMRI at insulating boundaries is attained when the magnetic field is nearly azimuthal, which produces a 
short-wave instability as in the WKB approximation. The short-wave character of the instability explains why the 
nonlocal solution obeys the Liu limit and why it may also apply to other gap widths rather than just the case of A = 2 
considered in this study. Moreover, the short wave length implies that the absolute instability, which formally applies 
to the systems of infinite axial extension, may be relevant also for relatively thin disks. This is because the absolute 
instability, as discussed at the beginning of Sec. IIVB1 represents the limit of the global instability with respect to the 
wave packets of finite length. The wave packets formed by two short waves with slightly different wave numbers can 
fit in the relatively thin disks. 

In the case of perfectly conducting boundaries, the instability can extend sufficiently far beyond the Rayleigh line 
to reach the astrophysically relevant Keplerian velocity profile. However, this is the case only for the convective 
instability mode, which is not self-sustained and thus requires an external excitation to be effective. Moreover, as the 
ratio of rotation rates increases further beyond the Rayleigh line, this instability mode shifts to ever longer waves. 
First, due to its long wave length this mode is not only inherently nonlocal but also unable to fit within limited 
thickness of accretion disks. Second, the critical interaction parameter for this mode approaching zero implies that 
its growth rate approaches zero too. As a result, this weak instability can be suppressed by a finite viscous damping, 
which also limits its extension beyond the Rayleigh line [Tol [2l|. Given these three basic constraints, the HMRI at 
perfectly conduction boundaries seems of a little astrophysical relevance, too. 

Note that the Liu limit appears only in the inductionless approximation, which, however, captures the essence of 
HMRI. At nonzero Pm, instability can extend beyond the Liu limit [8]. Obviously, the HMRI can extend similarly to 
the SMRI up to the Velikhov-Chandrasekhar limit corresponding to the solid-body rotation. However, this extension 
is due to the SMRI, in which HMRI merges with at a sufficiently high Reynolds number [2?]]. For small Pm, the 
merging point, at which the HMRI loses its identity by turning into the SMRI, is located only slightly above the Liu 
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limit [28]. Thus, the HMRI cannot affect the Keplerian velocity profiles even at finite Pm. 
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